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Abstract. Calculations of microscopic optical potentials (OP's) (their real and imaginary parts) are per- 
formed to analyze the 6 He+p elastic scattering data at a few tens of MeV/nucleon (MeV/N). The OP's 
and the cross sections are calculated using three model densities of 6 He. Effects of the regularization of 
the NN forces and their dependence on nuclear density are investigated. Also, the role of the spin-orbit 
terms and of the non-linearity in the calculations of the OP's, as well as effects of their renormalization 
are studied. The sensitivity of the cross sections to the nuclear densities was tested and one of them that 
gives a better agreement with the data was chosen. 
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1 Introduction 

The basic characteristics of the exotic nuclei, such as their 
charge and matter distributions have been tested, in par- 
ticular, by studying differential and total reaction cross 
sections of the proton scattering on exotic nuclei in inverse 
kinematics. Nowadays, a substantial amount of experi- 
mental data exists for the cross sections of 6 He+p elastic 
scattering at different energies. For example, proton elas- 
tic scattering angular distributions were measured at inci- 
dent energies less than 100 MeV/N for 6 He, namely 25.2 
[TPPPl- 38.3 [5], 41.6 [5117115] . and 71 MeV/N PIP] and 
at energy of 700 MeV/N for He and Li isotopes (e.g. refs. 
[TT1IT21IT51IT1IT5] V The analyses of differential and total re- 
action cross sections have been performed using different 
phcnomcnological and microscopic methods and models of 
nuclear structure (see, e.g. refs. ^^^j^^^MM^ 
[T7lfl8. 19, 20,21]). We note also the microscopic analysis 
(refs. [2l p 3 l l2 ^ l2l p 6 T l2"7]) based on the "coordinate-space 
g-matrix folding method" [28 (and some modifications 
[29 ), where non-local OP is obtained using a folding of a 
local medium-dependent NN effective interaction with the 
target ground-state mixed density. In some of the calcula- 
tions (e.g. [9, 10, 16, 17 ) the eikonal approach using proton 
and neutron density distributions, as well as parametrized 
NN total cross section have been used. It has been ac- 



a e-mail: gaidarov@inrne.bas.bg Present address: Insert the 
address here if needed 



cepted that for energies larger than 500 MeV/N the mul- 
tiple scattering diffraction theory developed in [301131] (the 
Glauber theory) is a relevant method to study charge and 
matter distributions from proton elastic scattering data 

The experimental information on cross sections of 
6 He+p elastic scattering requires for its adequate descrip- 
tion a development of the respective microscopic meth- 
ods of their analysis. These methods give an opportunity, 
first, to distinguish between different models for the ex- 
otic 6 He as a nucleus with a halo with two neutrons and, 
secondly, to test the attempts of the folding approach for 
constructions of optical potentials. The latter include the 
understanding of the role of various components of the 
OP and the necessity to introduce fitting parameters. A 
number of works has been devoted to calculations of OP's 
using the folding approach (see, e.g. [20,21,22,23], 34.35, 
36,37,38 ). For instance, the real parts of OP's for cal- 
culating the 6 He+p, 6 He+ 4 He (E lab =lbl MeV) [20] and 
6 He+p, 8 He+p (E mc < 100 MeV) [2T] elastic differential 
cross sections have been obtained microscopically using 
realistic M3Y-Paris effective NN interaction [35 , 38 , 39J to- 
gether with the Tanihata et al. proton and neutron densi- 
ties of the helium isotopes [40] in refs. [201121] and also with 
the densities of the cluster-orbital shell-model approxima- 
tion (COSMA) [5 ]fT0TIT8p ^ ] . In [21] a comparison of the 
obtained results has been performed to those from the 
alpha-core approach with the complex and fully non-local 
effective interaction [22] and also with the non-core model 
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based on the large-scale shell-model (LSSM) calculations 
(refs. [7]l3"6"l[37] and references therein). It was shown that 
the elastic scattering is a good tool to distinguish between 
different density distributions [21] . Usually, in the usage of 
the complex OP's for analyses of the differential cross sec- 
tions, their imaginary part and the spin-orbit terms have 
been determined in a phenomenological way, and then the 
OP's include a number of fitting parameters. The question 
to optimize this number in the analyses of the experimen- 
tal data is not usually considered. 

The main aim of the present work is to calculate dif- 
ferential cross sections of elastic 6 He+p scattering at dif- 
ferent energies studying the possibility to describe the ex- 
isting experimental data by using a minimal number of 
fitting parameters. We note that for this purpose we use 
the so-called high-energy approximation (HEA) OP [41] . 
Its form can be considered as a microscopic folding of the 
densities of the colliding nuclei NN scattering amplitude. 
There are no free fitting parameters in this OP and its 
dependence on the energy is included in the input data 
of the NN scattering amplitude and the total cross sec- 
tion. It is generally believed that the Glauber (eikonal) 
approximation is reasonable at energies of hundreds MeV 
and higher. However, beginning from the work [42] the 
method of HEA has been usually modified by replacing 
the eikonal straight-line trajectory at an impact parame- 
ter b by the parameter b c that corresponds to the distance 
of the closest approach of the projectile in the Coulomb 
plus the nuclear potential. For the medium and heavy nu- 
clei the Coulomb distortion dominates and successful ap- 
plications of the approach were demonstrated, firstly, in 
[43] and later on in many papers, e.g. in 41,44,45,46] at 
low (> 10 McV/N) and intermediate energies to the de- 
scription of the data on the differential elastic and total 
reaction cross sections of various projectile ions and target 
nuclei. The distortion effect caused by the nuclear poten- 
tial is important mainly for lightest nuclei as shown in 
[47] , In the last decade such a modified Glauber method 
turned out to be rather effective and was employed in 
many works devoted to analyses of nucleus-nucleus scat- 
tering processes. In the present paper, to avoid limitations 
in the modified HEA formulae, we account for distortion 
effects by computing cross sections using the DWUCK4 
code of numerical solving of the Schrodinger equation. At 
the same time we use (similarly to ref. [48]) the micro- 
scopic HEA imaginary part of the OP obtained in ref. [41] 
that yields the same eikonal phase as that given in the 
optical limit of the Glauber microscopic theory of mul- 
tiple scattering of complex systems. So, one of the main 
aims of our work is to establish the limits of the applica- 
bility of the HEA OP for calculations of differential elastic 
cross sections of the 6 He+p scattering for different regions 
of angles and incident energies. Together with the HEA 
imaginary potential we tested the OP whose form coin- 
cides with that of the real part of the standard folding OP 
[34,35,38]. The latter includes an exchange term and, cor- 
respondingly, the non-linear effects in calculations of the 
potential. In the calculations we pay attention to the role 
of the importance of the contribution of various physical 



quantities and features, such as microscopically obtained 
spin-orbit forces and regularization of the NN forces used 
in folding calculations. Also, we consider effects of renor- 
malization of the real and imaginary parts of microscopic 
optical potentials, the differences or similarities of various 
models for the 6 He structure which are used in the de- 
scription of the experimental data on the cross section of 
the 6 He+p elastic scattering. 

The theoretical scheme for microscopical calculations 
of the real part of the OP's and cross sections is given 
in sect. 2. This sect, includes also some methodical calcu- 
lations. Section 3 is devoted to the OP within the HEA. 
The results of the calculations and the discussion are pre- 
sented in sect. 4. Section 5 includes the conclusions from 
the work. 



2 Basic relationships for calculations of the 
real part of the nucleon-nucleus optical 
potential 

2.1 Direct part of the real OP (Re OP) 

The real part of the nucleon-nucleus OP is assumed to be 
a result of a single folding of the effective NN potential 
with the nuclear density, i.e. this is a particular case of 
the double-folding [34 in which a S(ri) function has to be 
used for the density of the incoming particle p(ri). Then 
the direct part of the Re OP (V D ) has the following form 
of the isoscalar (IS)- and isovector (IV)- contributions, 
correspondingly: 



VP s ir) = I p 2 {v 2 )g{E)F{p 2 )v» Q {s)d A r 2 , 



D i 



(1) 



Vj D v (r) 



where s = r + r 2 , 



6p 2 (r 2 )g(E)F(p 2 )v»(s)d 6 r 2 , (2) 



p2(r 2 ) = p 2 ,p(r 2 , p ) + p 2 , n (r 2 , n ), (3) 

Sp 2 {r 2 ) = p 2 , P (r 2 , p ) - P2,n(r2,„). (4) 

Here p 2lP (r 2 . p ) and p 2t n(Y 2 , n ) are the proton and neu- 
tron densities in the target nucleus. In eqs. ([1]) and ([2]) 
g(E) = 1 — 0.003E 1 represents the energy dependence of 
the effective NN interaction while F{p 2 ) contains its den- 
sity dependence. Following ref. [38] we use its form for the 
CDM3Y6 effective Paris potential: 



F(p) =C 1 + ae 



7p(r) 



(5) 



where C = 0.2658, a = 3.8033, = 1.4099 fm 3 , 7 = 4.0 
fm 3 . 

The effective NN interaction t^( i) in eqs. ([1]) and (J2| 
has included the isoscalar and isovector components of the 
direct part of the M3Y interaction based on the results 
of the g-matrix calculations using the Paris NN potential 
[551155] . The M3Y potentials which are used in the fold- 
ing calculations of OP's are sums of Yukawa- type terms 
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exp(— fir) /(fir). Using eqs. (U)-© one can obtain the fol- 
lowing forms of the direct part of the isoscalar Re OP 
expressed by integrals in the coordinate and momentum 
space, correspondingly: 



(6) 



Vps( r ) = Cg(E) J |p 2 (r 2 ) + ag 2 (v 2 ) - 7P2(r 2 ) 
x v£ (s)d 3 r 2l 



VP s (r) 



Cg(E) 
2tt 2 



p 2 (q) + ag 2 (q) ~ 192(g) 



x v oo(q) Jo(qr)q dq, 



(7) 



where Q 2 (r 2 ), g 2 (r 2 ) and their Fourier transform have the 
forms: 

Q2(r 2 ) = P2 (r 2 )e-^^\ (8) 
£ 2 (r 2 ) = h(r 2 )] 2 , (9) 



p(q)= J e lclr p(r)d 3 r = Air p(r)j (qr)i 



2 dr. 



(10) 



Similarly, exchanging p 2 by 8p 2 [eq. one can obtain 
the isovector part Vjy of the direct part Re OP. 

2.2 Exchange part of the real OP (Re OP) 

In contrast to the case of the double-folding potential 
where integration over the coordinates (ri) of the nucle- 
ons in the incoming nucleus takes place, in the case of the 
nucleon-nucleus interaction one can obtain: 

V I E s x (r)=g(E) J p 2 (r 2 , r 2 - s)F (p 2 (r 2 - s/2)) 

X vg x (s)j (k(r)s)d 3 r 2 . (11) 

For the density matrix p 2 (r 2 ,r 2 — s) in cq. (TTTj) we use 
the approximation for the calculation of the knock-on ex- 
change term of the folded potential from |49j which pre- 
serves the first term of the expansion given in [50j : 

p 2 (r 2 ,r 2 -s) ~ p 2 (|r 2 -s/2|)j 1 (fc Fi2 (|r 2 -s/2|)- S ). (12) 
In eqs. (dJ) and (fT2"]) : 



ji(x) = —ji(x) = — (smx — xcosx), 



(13) 



d 3 s = s 2 ds sm8d6dcf), with 9, cj> being angles which de- 
termine the vector s/2 with respect to the fixed vector r. 
Integrating over dip gives the factor 2tt. We separate the 
integral over dx = sin6>d# = — dcos#: 



G l S (r,s) = j p 2 (x 2 (r,s,x)) ji(k F}2 (x 2 (r,s,x))) 
x F(p 2 (x 2 (r,s,x))) dx, 



where 

x 2 (r,s,x) = |r-s/2| 



1/2 



(17) 



(18) 



Finally, the isoscalar part of the exchange Re OP of the 
nucleon-nucleus interaction has the form: 



Vt s X (r) =2ng(E) 



x j G (r, s)v* x (s)j (k(r)s) s'ds. (19) 

In eqs. (TTTj) and (|19p . Vq Q x (s) is the isoscalar part of the 
exchange contribution to the effective NN interaction and 
the local momentum of the incident nucleon in the field of 
the Coulomb and nuclear potential (Re OP) is [5T] : 



k 2 = (2m/h 2 ) [E c . m . - Vc(r) - V(r)} 



l + A 2 



-4, 



(20) 



In eqs. (fT2"]) , P^|) and ([TTJ) kp, 2 defines the average relative 
momentum [49ll51j : 



'-F.2 



(r) 



5 



r(p) ~ ~^ 2 p(r) 



1/2 



(21) 



where we choose further the extended Thomas-Fermi ap- 
proximation 52,52] for the kinetic energy density 



r(p) 



3 



r q (p q )^^(^ 2 Y^[p q (r)] 

\Vp q (r)\ 2 | V 2 Pq (r) ^ 
36p q (r) 3 



5/3 



(22) 



valid for each kind of particles q = n,p. 

The isovector part of the exchange Re OP can be ob- 
tained from eqs. (fTF]) and (fT9|) exchanging p 2 by Sp 2 . 



F(p 2 ) = C [l + ae-l 3 ^ 2 -^ - 7P2 (r 2 - s/2)] . (14) 

In our case r+r 2 = s, therefore |r 2 — s/2| = |r— s/2| = |x 2 | 
and making the substitution d 3 r 2 =d 3 s at fixed r, eq. (jlip 
can be rewritten in the form: 



V I E s x (r)=g(E) / h 2 (r - s/2, s)F (p 2 (\r - s/2|)) 



where 



x v& x (s)j (k(r)s)d 6 s, 



h 2 (x 2 ,s) = p 2 (x 2 )ji(/cF, 2 (x 2 ) • s) 



(15) 
(16) 



2.3 "Regularization" of the effective NISI interaction in 
the form of Yukawa potential 

An important aspect of using the M3Y effective NN po- 
tentials in their traditional form 



v(s) 



E 

3 



N 



exp(-fij\s\ 

flj\s\ 



(23) 



is related to their regularization due to the singularities. 
The latter aims to exclude the singularities at the point 
|s| = 0, which have no physical meaning. Formally, this 
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singularity might affect the consideration of the nucleon- 
nuclcon scattering, but it does not cause any difficulties 
in the calculations of the folding integral [e.g. eqs. (|6|) 
and (fl~5]) ] to generate the nucleon-nucleus and/or nucleus- 
nucleus potential. Nevertheless, "regularized" M3Y poten- 
tials of the NN interaction which do not contain the men- 
tioned singularity have been used in some works {e.g. [21[ 
I54j ) and in our work we would like to consider this point 
in more details. Note that the physical reason for that 
kind of problems is related to the breaking of the meson 
theory at very short range due to the extended structure 
of nucleons [55] . That is why the one-boson-exchange po- 
tentials are usually regularized by introducing, e.g., the 
monopole-, dipole-, exponential cut-off form factors [56] . 
as well as that suggested in ref . [57] . As shown in [58] , the 
most important feature of these form factors is that they 
have the same principal range, and that the overall results 
are insensitive to their detailed shape. In refs. [211154] cut- 
off form factor was taken to be equal to the proton form 
factor p p (q), parametrized according to the experimental 
data in |59] as a sum of Gaussian functions. As a test 
of the role of the regularization, in the present work we 
take the same expression, and then use the corresponding 
"smoothing" function in a coordinate space, which has a 
meaning of the density distribution of the incoming proton 



P P ( r ) =X/ 



1 



> r 2)3/2 



exp 



(24) 



where = 1, a 1= 0.506373, a 2 =0.327922, a 3 =0. 165705, 
r\ = 0.431566 fm 2 , r\ = 0.139140 fm 2 , r\ = 1.525540 fm 2 , 
and (r 2 )=0. 77542 fm 2 . Then, the regularized NN potential 
is determined by 



(si) = / p p (r 1 )v(s)d n, 



si = ri 



(25) 



Rewriting this integral in the momentum representation 
one gets the expression used in [2"TU5"4] instead of v(s) 



Vreg(si) 



1 



B "^'~ (2tt) 3 
where the proton form factor 



p p (q)v(q)e l ^d 3 q, (26) 



2 2 

q rf 



PM) = J e- iqri p p (n)d 3 n = exp (^-^ ) (27) 



and 



„r 47T 1 



(28) 



On the other hand, substituting the "regularized" v reg (si) 
potential ([2"S"]) in eq. (1) instead of v(s) leads to the ex- 
pression 

VP s {r) = J p p (r 1 )g(E)F(p)vS ) ( S )p 2 (r 2 )d 3 r 1 d 3 r 2 , 

ri+s = r + r 2 , (29) 

whose form coincides with the definition of the double- 
folding potential of interaction of two colliding complex 
systems having densities p p and p 2 . 



2.4 Methodical calculations 

In this subsection we will present results of our calcula- 
tions of the 6 He+p elastic scattering differential cross sec- 
tions studying the role of various factors, such as: i) the 
choice of the density distribution of 6 He, ii) effects of the 
regularization, and iii) the effect of the spin-orbit term in 
its dependence on the interaction potential. 

In the calculations we use the following three density 
distributions of 6 He: 

i) the point-nucleon density 



y 

Ppoint 



+ 



7T 3 / 2 

1 X 



1 

-r exp 

cr 



6 3 3 



exp 



(30) 



applied by Tanihata et al. [60] for a comparison of the 
measured total reaction cross section of 6 Hc+ 12 C at 800 
A MeV with the respective expression from [61] obtained 
there in the optical limit of the Glauber theory. In (|3H]l 
X = Z.N, and the parameter values of a and b are deter- 
mined from 



2 *2 I i 

a = a I 1 



1 

A) ' 



6 2 = b* 2 { 1 ' 



A ' 



(31) 



where a* — 1.53 fm, b* = 2.24 fm, and hence a = 1.40 
fm, b = 2.04 fm. Thus, the rms radius of the point-proton 
density of 6 He is equal to 1.72 fm. 

ii) the COSMA point-nucleon density [TH] which has 
the same analytical form as eq. (130[) . but with the param- 
eter values a = 1.55 fm and b = 2.24 fm [9][T0]. and the 
rms radius of the point-proton density is equal to 1.89 fm. 

We should emphasize that both Tanihata and COSMA 
densities have a Gaussian asymptotic behavior which is 
not a realistic one at high q. That is why we consider also 
other more realistic proton and neutron densities of 6 He: 

iii) the LSSM densities obtained in a complete 4hu> 
shell-model space [36] using Woods-Saxon (WS) single- 
particle wave function basis with realistic exponential 
asymptotic behavior. 

In fig. [T] we present the total density distribution of 
6 Hc: p(r) = p p (r)+p n (r) (in logarithmic and linear scale), 
as well as the point-proton and point-neutron density dis- 
tributions of Tanihata et al, COSMA and LSSM. One can 
see that the Tanihata and COSMA densities have a Gaus- 
sian slope while the LSSM tail occurs too higher and goes 
for the neutrons to larger values of r than for protons. As 
known, the differences between the densities for smaller 
values of r (r < R, R being the linear size of the nucleus) 
can be revealed mostly in the nucleon-nucleus scattering, 
but in heavy ion collisions one can study better the asymp- 
totic region. 

The role of the regularization (described in subsect. 
2.3) is shown in fig. H]by a comparison of the 6 He+p elas- 
tic scattering differential cross sections at different ener- 
gies: 25.2, 41.6, and 71 MeV/N. As can be seen, generally 
the effect of the regularization is rather weak, but it in- 
creases when the energy and the scattering angle increase. 
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lt starts to be seen at smaller angles when the energy in- 
creases. Generally, checking the role of the regularization 
we conclude that it is not necessary to include it in the 
cases of folding OP's in the nucleon- nucleus scattering. 



10° 



E=25.2 MeV/N 




e c. m . [ de 9l 

Fig. 2. Elastic He+p scattering cross sections at different 
energies computed by using the microscopically folded real 
OP and the imaginary part taken in the same form (V F = 
V D + V EX ) calculated using the LSSM density of 6 He with 
(dashed line) and without (solid line) regularization of the M3Y 
effective NN interaction. 



In fig. [3] we show the role of the spin-orbit term in the 
cross section of elastic 6 He+p scattering at a given energy 
E = 41.6 MeV/N. The spin-orbit term is taken to be in 
the form: 

»~* (;) < 32 > 



where f(r) is the form of the Re OP, and fm 2 . The 

calculations are given for three cases: i) when N so — 0.5 
and f(r) is exchanged by the microscopic Re OV=V mlcro (r) 
in MeV; ii) N so = 6.2 MeV and f(r) is the WS form fac- 
tor taken from [62] . and iii) without spin-orbit term, i.e. 
when N so = 0. 

One can see that the cross sections in the cases i) and 
ii) are very close to each other. Our analysis showed that, 
generally, there is not a strong dependence on the shape of 
/(r). From the other side, however, the effect of the spin- 
orbit term in the considered case is important at larger 
angles (> 60°). 



JO 



a 
To 




100 



Fig. 3. Elastic He+p scattering cross section at energy E — 
41.6 MeV/N calculated using the microscopic Re OP (V F ) in 
the spin-orbit term (l/r)dV^ F /dr (solid line) and also using the 
WS Re OP in 12.4(l/r)d/ ws /dr (dotted line). The dashed line 
shows calculations without spin-orbit term. The LSSM density 
is taken for 6 He and the Im OP has the form of V F . 



3 Optical potential within the high-energy 
approximation 

As known, the real part of the optical potential can be cal- 
culated microscopically in the standard form of the single- 
or double-folded integral. In principle, the physical nature 
of the real and imaginary parts of OP is different. The 
imaginary part is related to the flux loss at the transition 
of the particles from the elastic to the inelastic and reac- 
tion channels which depends on both the structure of the 
colliding nuclei and reaction mechanisms. This makes it 
difficult to construct practically a convenient theory of the 
Im OP especially for nucleus- nucleus scattering (see, e.g., 



6 



K.V. Lukyanov et al.\ Calculations of 6 He+p elastic scattering cross sections 



[63 and refs. therein). However, in the case of more simple 
proton-nucleus scattering one can get fairly good applica- 
tions (see, e.g., [64]) when using a single-folding pseudo- 
potential [34] multiplied by the fitted complex renormal- 
ization factor (Nr + i/Vj) to obtain the complex potential 
with the unit shape of the real and imaginary parts. On 
the other hand, in the case of heavy ion scattering, many 
applications were made when the real part of OP was mi- 
croscopically calculated while an imaginary part was taken 
in the WS form with three or more fitted parameters. 

In this paper we intend to test the so-called HEA op- 
tical potential (at least its imaginary part) to explain the 
available data on the 6 He+p differential cross sections at 
a few tens of MeV/N. As already mentioned in the Intro- 
duction, the generalization of the eikonal method made it 
possible to use it at relatively low energies [41,42,43,44. 
45,46,47 . In our work, on the base of the eikonal phases 
we reconstruct the equivalent potential including its imag- 
inary part. The calculations of the cross sections are per- 
formed by numerical integration of the Schrddinger equa- 
tion by means of the code DWUCK4 using all interactions 
obtained (Coulomb plus nuclear OP), but not using the 
HEA scheme for cross section calculations. We note that 
using this method we aim to establish the limits of the 
applicability of the imaginary part of the HEA optical po- 
tential. We also note that, firstly, in the calculations of the 
HEA potential, similarly to ref. [4"TH4"5] . we do not neglect 
the parallel contribution of the transferred momentum (in 
contrast to the Glauber method). As a result, the integral 
in the expression for the potential includes spherical Bessel 
functions, instead of the cylindrical ones (which appear in 
the expressions for the eikonal HEA). Secondly, this makes 
softer the limitation of the applicability of the method. We 
mention that in this case the region of the applicability 
of the standard HEA for E 3> |f^(-R)| and small angles is 
transformed to E < \U{R)\ and <d < [y/2/kR+\U(R)\/E], 
with R being the radius of the potential at its periphery 
(e.g. at the half-depth of the nuclear potential), where 
the absorption is quite strong [65j . Previously, the micro- 
scopic HEA OP was successfully applied in refs. [41, 66, 67J 
to describe the 16,17 0+A elastic scattering data at about 
hundred of MeV/N from [68,46]. Here we would like to 
give an estimation for the angles of the application of our 
approach. For the 6 He+p scattering at i?; a b=40 MeV/N 
and R ss 2 fm, it follows from the expression given above 
that i? < 50°, i.e. this is the applicability region of angles 
considered in our calculations. The basic component of 
HEA is the eikonal phase <P(b) which depends on the im- 
pact parameter of the collision b. The amplitude and the 
cross sections of scattering and reactions are expressed 
by means of ^(6). In the phenomenological approach the 
phase is given by the integral along the axis z of the inter- 
action potential. On the other side, this phase was derived 
in the microscopic level for the proton-nucleus scattering 
in the Glauber theory [3"Dll3"T] and generalized later to the 
nucleus-nucleus scattering in 69 , 70J . In the so-called "op- 
tical limit" of this theory an explicit form of the micro- 
scopical phase <P(b) is expressed by the densities of collid- 
ing nuclei and of the amplitude of the NN scattering. By 



a comparison of both the eikonal and microscopic expres- 
sions for the phase it became possible in [UJ and [JS] to 
obtain the explicit form of the OP in HEA, which gives a 
description of the nucleus- nucleus scattering, being equiv- 
alent to that from the microscopical approach: 



iW 



H 



(33) 



where 



V H (r) 



hv 



T°~NNC*NN 



(2tt)2 

Pi(q)P2(q)fN(q)jo(qr)q 2 dq, (34) 



W H (r) = 



hv 



T&NN 



(27T) 2 

oo 

Pi(q)P2(q)fN(q)jo(qr)q 2 dq, (35) 



with v being the velocity of the nucleus- nucleus relative 
motion, p\i2)(q) being the form factors corresponding to 
the point-like nuclcon density distributions of the nuclei 
and /at (q) being the amplitude of the NN scattering which 
depends on the transfer momentum q (see, e.g. [71]). The 
quantities cr/v n and a nn are the averaged over the isospins 
total NN cross section and the ratio of the real to imagi- 
nary part of the scattering amplitude at zero angle of free 
nucleons. They can be obtained from the data on the mu- 
tual scattering of nucleons. For instance, the parametriza- 
tions of the energy dependence of unn |44j and o.nn in 
the interval e ia b = 10 MeV-r- 1 GeV [JJ] are known: 



0~NN 



Z\Z 2 a vv + NiN 2 cr nn + Qcr n 



AiA 2 



(36) 



C = Z X N 2 + N X Z 2 



a np = (-70.67 - 18.18/T 1 + 25.26/T 2 + 113.85/3) 
x fm{np), (37) 



= (13.73- 15. 04/T 1 
x f m (nn), 



8.76/3- 2 + 68.67/3 4 



where 



931.5 



Slab + 931.5 



(38) 



(39) 



is the ratio of the relative to the light velocities, e; a (, = 
E/Ai is the energy (in MeV) of a nucleon in the inci- 
dent nucleus in the laboratory system Q, and cross sec- 
tions are given in mb. The factors f m (np) and f m (nn) 
are introduced to correct the dependence of cross sections 
on the energy and on the density of nuclear matter. The 
in-medium NN interaction has been widely investigated. 
For instance, in [73] numerical calculations of the total 
NN cross sections have been performed on the basis of 



In our case one has Z\ = Ai = 1, Ni = 0. 
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the Dirac-Brueckner theory of nuclear matter and their 
parametrization has been given in ref. [74j to obtain the 
following correcting factors: 



fm{np) = 



f m (nn) 



1 



20.88e° a °V° 2 



1 + 35.86p 



1.90 



1 



7.772 £ ° Q °y- 48 



1 



(40) 



(41) 



3 . In the case 
The increase 



18.01/9 1 - 46 

In eqs. (j4"D|) and (|4"Tj) the densities are in fm 
of free nucleons (p=0) f m {np)=f m {nn)=l 
of the density leads to decrease of these factors and, cor- 
respondingly, of the cross sections. 

Here we present also the ratio of the real to imaginary 
part of the NN amplitude at zero angle averaged over the 
nuclear isospins and parametrized in [72 : 



OtNN = 



Z\Z 2 a P p(j pp + NxN 2 a 



np u np 



Z 1 Z 2 a pv + NiN 2 (J ri 



(42) 



= 0.0078 + 0.1762 
-0.0301 + 0.2148 




0.01436e, ob , 
0.055L 



(43) 

*np - - U .uuu.L -r v.^^ySlab - U.UDDi£ iafc . (44) 

The results of calculations of the 6 He+p elastic scat- 
tering cross sections at energy 41.6 MeV/N using the OP 
in the form (|33p are presented in fig. 0] with and without 
spin-orbita term of the form of 0.5A 2 .(d/dr)I/ mlcro . It is 
seen a better agreement with the data up to i9 ~ 40° in 
the case when the spin-orbit term is included. The LSSM 
densities of 6 He have been used in the calculations. Here 
we would like to mention the question whether the data 
can be reasonably fitted either by introducing a spin-orbit 
force, or by renormalizing the potential. For the folding 
potential these two ways are not equivalent because of the 
role of the spin-orbit term at large angles. However, in 
the case of the HEA, responsible for the area of compa- 
rably small angles, it seems likely that one could fit the 
data using either a spin-orbit term or a renormalization. 
We note, however, that our main goal is to show that the 
HEA OP's can be used in their domain of validity with 
no additional free parameters; this is certainly true with 
the inclusion of the spin-orbit term at angles smaller than 
40°. 

In the next sect., following the basic theoretical scheme 
given in sects. 2 and 3, we present the results of our cal- 
culations. For calculations of Re OP we use the effective 
M3Y interaction based on the Paris NN potential with a 
density dependence in the form CDM3Y6 from [35] with- 
out "regularization" . Also we account for the spin-orbit 
term ~dy mlcro /dr and use the LSSM proton and neutron 
densities of 6 He. Besides, we include in the calculations the 
imaginary part of the OP's obtained within the HEA. 



4 Results of calculations and discussion 

In this sect, we present our calculations of the cross sec- 
tions of 6 Hc+p elastic scattering at different energies aim- 
ing to study: i) effects of the different behavior of the 



-S2 



Ci 

"D 
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U =V n +iW n 
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He+p 
E=41.6 MeV/N 



e c. m . [ de 9l 

Fig. 4. Elastic 6 He+p scattering cross section at energy E — 
41.6 MeV/N calculated by using U op t = V H + iW H with (solid 
line) and without (dashed line) spin-orbit term (l/r)dV H /dr. 
The 6 He LSSM density is used. Experimental data are from 
refs. [6l|8]. 



model densities of 6 Hc; ii) the possibility to use the HEA 
optical potentials for different energies and angles and 
also the microscopic OP with the folded real part V F = 
yD + yEX anc j tne jjea imaginary part W H ; hi) the role 
of the renormalization of the depths of the real and imag- 
inary parts of OP's, and iv) the in-medium effects of the 
effective NN interaction on the microscopically calculated 
Re OP. 

We start this sect, with calculations of the 6 He+p elas- 
tic cross sections using the real part of OP calculated 
within the folding approach (V F ), as well as the real (V H ) 
and imaginary (W H ) parts of the HEA OP's. We intro- 
duce two renormalization parameters Nr and Nj (already 
mentioned above) and consider the following three types 
for the OP's: 



(A) 
(B) 



Kt = N n vH 



iNfW H , 



K = NiV F + iNf W», 

K = N R yF + iN Fv F . 



(45) 

(46) 
(47) 



As can be seen, in case (A) we use both real and imagi- 
nary parts from the HEA calculations of OP's; in case (B) 
we take V F , the folded real part of the microscopic OP 
where the exchange term is included, and the imaginary 
part is applied in the form (|35|) of the HEA OP, while in 
case (C) we use the microscopically folded real and imagi- 
nary parts in the form of V F , and thus they have the same 
shape. 

In fig. [5] are presented the results of calculations of the 
6 He+p elastic cross sections (for energy E=A1.6 MeV/N) 
with the fixed value N/=l for all three cases (A), (B) 
and (C). The comparison with the data is performed for 
two values of Nr for each case, namely 0.53 and 1.00 for 
case (A), 0.85 and 1.00 for cases (B) and (C). One can 
see a good agreement with the data using Nr=1.00 for 
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the cases (A) and (B) and using Nr = 0.85 for case (C). 
It can be concluded that for both cases (A) and (B) the 
renormalization is not necessary. As in the case shown in 
fig. HI we note that the agreement obtained by using the 
HEA is for angles smaller than 40°. 
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Fig. 5. Elastic 6 He+p scattering cross section at energy E — 
41.6 MeV/N calculated using different OP's (cases (A), (B) 
and (C) from the text) for various values of the renormalization 
parameters Nr (Nj = 1). The LSSM density of 6 He is applied. 
Experimental data are from refs. [6]|8]. 



The latter concerns also the results of the calcula- 
tions of 6 Hc+p elastic scattering (£7=41.6 MeV/N) in case 
(B) using three densities of 6 He given by Tanihata et ai, 
COSMA and from the LSSM given in fig.[H It can be seen 
a good agreement with the empirical data up to d ~ 35°. 
Deviations of the results with the use of the COSMA den- 
sity from the other two cases start at angles larger than 
40°. In these calculations N R = Nj = 1. 

In fig. [7]we present the real (V F ) and imaginary (W H ) 
parts of the OP's (case (B)) calculated using the LSSM 
density of 6 He for three different energies: 25.2, 41.6, and 
71 MeV/N. They are given without renormalization (i.e. 
Nfi=Ni—l) and are illustrated in logarithmic and linear 
scales. One can see the decrease of the potential depths 
with the increase of the energy. Here we would like to note 
that, as can be seen from fig. [71 at R ps 2 fm the values of 
the potentials for E iab = 40 MeV/N are U{R) ps 15 - 20 
MeV, and thus the limit condition for the HEA E > U(R) 
is fulfilled. This explains the applicability of the HEA at 
such values of the energy. We should mention the larger 
depth of the HEA imaginary parts (W H ) of the OP for the 
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Fig. 6. Elastic 6 He+p scattering cross section at energy E — 
41.6 MeV/N calculated using U op t = V F + iW H and Tanihata 
(dashed line), COSMA (dotted line) and LSSM (solid line) 
densities of 6 He. Experimental data are from [5][S]. 




r [fm] r [fm] 

Fig. 7. Microscopic real part (V F ) of OP ((a) and (a')) and 
HEA imaginary part (W H ) ((b) and (b')) calculated using the 
LSSM density of 6 He for energies E = 25.2 (solid lines), 41.6 
(dashed lines) and 71 MeV/N (dotted lines). 



case of 25.2 MeV/N seen in fig.[7J As we will see below, the 
result of the HEA for W H in this case affects significantly 
the cross section. 

The three densities are used to calculate the cross sec- 
tions of 6 He+p elastic cross sections for three energies 
(25.2, 41.6, and 71 MeV/N) shown in fig. M One can see 
the fairly good agreement with the experimental data of 
the results with the LSSM density for 6 He for energies 
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corresponding depths NnVp T ~°' > and NjW^~ v> of the real 
and imaginary parts of OP's. The values of Nr and TV/ 
were chosen starting from the values Nr = 1 and Nj = 1 
and decreasing them gradually in order to achieve a rea- 
sonable fit to the experimental data. In our opinion, the 
obtained values of Nr and Nj still do not reveal some 
regular change with the increase of the energy. 

In fig.[5]a particular attention is paid to the case of the 
elastic 6 He+p cross section at energy of 25.2 MeV/N cal- 
culated using the LSSM density for the 6 He nucleus. One 
can see that, in contrast to the case of larger energies, in 
this case smaller values of Nr (0.35) and especially of Nj 
(0.03) are necessary for a better agreement with the ex- 
perimental data. This concerns the slope of the cross sec- 
tion for angles i9 c . m . between 70° and 120°. We note that 
the results shown in fig. [9] for the energy 25.2 MeV/N 
are for angles up to d c . m . ss 120°. The necessity to use 
much smaller value of the renormalization parameter Ni 
(Ni = 0.03) to fit the data for large angles is related to 
the large value of the depth of the imaginary part of the 
OP obtained within the HEA (W H ) for this al- 
ready mentioned above. This already shows the limitation 
of the approach (the case (B), Eq. ([46]) ) for small ener- 
gies (< 25 MeV/N) and large angle values. In this case 
the condition E > \U(R)\ is not already fulfilled. Along 
this line we should mention the Ref. [75] in which a micro- 
scopical pseudo-folding potential (without accounting for 
the exchange) V rnicro = (Nr + iN^Vfoiding was used for 
calculations of p+ 18 Ne and p+ 18 scattering at energies 
24.5 and 30 MeV/N. It was shown a good agreement with 
the data with TV/ = for 18 Ne and 7V 7 = 0.006 for 18 0. 
These results confirm the necessity to use microscopical 
rather than phenomenological OP's also for low energy 
scattering. 

In fig. [10] we present the in-medium effect of the NN in- 
teraction on the clastic scattering cross sections of 6 Hc+p. 
The calculations were performed using the function F(p) 
in the form ([5]) and with F(p) = l (the case without in- 
medium effect). The effective M3Y interaction based on 
the Paris NN potential and the LSSM density of 6 Hc were 
used in the calculations. One can see the existence of small 
differences between the results using density dependent 
and independent effective NN interaction. The differences 
for the larger angles increase with the energy increase. 

Also, we studied the non-linear effects in calculations 
of the exchange part of OP. As can be seen, the factor 
jo(k(r) ■ s) takes place in the expressions for the V EX 
(see eqs. ([TT]) . ([TS"]) and ([T9]) V where the local momentum 
of the relative motion k(r) [eq. ([20]) ] is expressed by the 
OP. The calculations need an iteration procedure and this 
makes the task complicated. However, if k(r) ~ 0, then 
jo — 1 and this simplifies calculations. In fig.[TTJwe present 
the results for the 6 He+p elastic scattering cross section 
at different energies with and without accounting for the 
factor jo(fc(r) • s) in the above mentioned equations. As 
can be seen, the non-linearity effect is small for the energy 
E = 25.2 MeV/N, but it increases with the energy increase 
up to an order of magnitude for i3 > 40° for the energy E = 
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Fig. 8. Elastic 6 He+p scattering cross sections at different 
energies calculated using U opt = N R V F + iNiW" for various 
values of the renormalization parameters Nr and Ni giving a 
reasonable agreement with the data (presented in table[TJ). The 
used densities of 6 He are LSSM (solid line), Tanihata (dashed 
line) and COSMA (dotted line). Experimental data are taken 
for 25.2 [Tpp]. 41.6 0(8] and 71 MeV/N 0[TD]. 



41.6 and 71 MeV/N, in contrast to the results obtained 
with the other two densities for the energies 25.2 and 71 
MeV/N. In table[JJwe list the values of the renormalization 
parameters Nr and Nj that give a reasonable agreement 
with the data for the three energies and the three differ- 
ent densities of 6 He shown in fig. [S] One can see also the 
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Table 1. The optimal values of the renormalization parameters N R and Ni obtained by fitting the experimental data for 
the elastic 6 He+p cross sections. In the calculations U op t = N R V F + iNiW H and LSSM, Tanihata and COSMA densities for 
energies _B=25.2, 41.6, and 71 MeV/N are used (the results are shown in fig. [8}. The depths of the corresponding potentials (in 
MeV) are presented, as well. 



Energy 


25.2 


25.2 


41.6 


41.6 


71 


71 


Density 


N R 


Ni 


N R 


Ni 


N R 


Ni 


LSSM 


0.6 


0.8 


1.0 


1.0 


0.6 


1.0 


Tanihata 


1.0 


0.6 


1.0 


1.0 


1.0 


0.5 


COSMA 


1.0 


0.6 
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1.0 
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Fig. 9. Elastic 6 He+p scattering cross sections for E — 25.2 
MeV/N calculated by using the LSSM density for 6 He. The 
curves exhibit results for U op t = N R V F + iNiW H with dif- 
ferent values of N R (0.6-dashed, 0.5-dotted, 0.4-dash-dotted) 
and fixed value of AT/=0.1. The solid curve is for N R —0.35 and 
iV/=0.03. The experimental data l 2 3 are also given. 



71 MeV/N. These results show the necessity to perform 
calculations accounting for the non-linearity of the task. 



5 Conclusions 

The results of the present work can be summarized as 
follows: 

i) The optical potentials and cross sections of 6 He+p 
elastic scattering were calculated at three different ener- 
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Fig. 10. In-medium effect of the M3Y NN interaction on cal- 
culations of elastic 6 He+p scattering cross sections for dif- 
ferent energies. The Re OP (V F ) is calculated with F(p) = 
C(l + a cxp(— ftp) — 7p) (solid lines) and with F(p) = 1 (dashed 
lines). The imaginary part is calculated within the HEA (W H ). 



gies E = 25.2, 41.6, and 71 MeV/N. The following com- 
ponents of the OP's were used: 

a) the real part of the OP (V F ) calculated microscopi- 
cally using the folding procedure and M3Y effective inter- 
action based on the Paris NN potential; 

b) the real (V H ) and imaginary (W H ) parts of the OP 
calculated within the high-energy approximation (HEA); 

c) three different combinations of V F , V H and W H 
(cases (A), (B) and (C), eqs. (JSJ-fllZ])) were used for the 
OP U op t in calculations of the elastic 6 He+p cross sec- 
tions. The renormalization parameters Nr and Nj have 
been introduced and their role has been studied. The cross 
sections were calculated by numerical integration of the 
Schrodinger equation by means of DWUCK4 code using 
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Fig. 11. The non-linearity effect of elastic 6 He+p scattering 
cross section for different energies. The solid lines are the cal- 
culation results when Re OP (V F ) includes jo(k(r) ■ s) term in 
eq. (11) and the dashed lines are without this term. The LSSM 
density of 6 He is used. The imaginary part of OP is calculated 
within the HEA (W H ). 



all interactions obtained (Coulomb plus nuclear optical 
potential); 

d) three different model densities of protons and neu- 
trons in 6 He were used in the calculations: the phenomeno- 
logical ones in the form (|30p . parametrized by Tanihata et 
al., the same form (|30|) with parameters from the COSMA, 
and also the microscopically calculated density within the 
LSSM. 

ii) The results of our calculations show that the LSSM 
density of 6 He is the most preferable one because it leads 
to a better agreement with the data for the 6 He+p elastic 
scattering at the three energies. The physical reason for 
the latter is that the LSSM densities have more diffuse 
tails at larger r than the densities based on Gaussians. 
With the LSSM density and with the optical potentials of 
the form (B), namely U opt = N R V F + iN^" , we tried 
to choose the parameters Nr and Ni starting from the 
values Nr = 1 and JVj = 1 and decreasing them gradually 
in order to achieve a reasonable agreement of the calcu- 
lated cross sections with the available data. The obtained 
set of these parameters is the following one: Nr = 0.6, 
1.0 and 0.6 and Nj = 0.8, 1.0 and 1.0 for energies 25.2, 
41.6, and 71 MeV/N, respectively. We note that the use 
of the microscopic folding real part V F and of the HEA 
imaginary part W H leads to a good agreement with the 
data for 41.6 and 71 MeV/N, while the data at lowest 
energy 25.2 MeV/N are explained only on the qualitative 
level. As shown by the estimations presented, this is re- 
lated to the limitations of using the imaginary part of the 
HEA OP for energy smaller than around 25 MeV/N due 
to the fact that the potentials do not fulfill the applicabil- 
ity condition E > \U(R)\. In this case the large value of 
the depth of the Im OP obtained in the HEA W H has to 
be strongly reduced {e.g. using in our case Nj = 0.03) in 
order to achieve a reasonable agreement with the data for 
the energy 25.2 MeV/N. 



iii) It was shown that the effect of the regularization 
of the M3Y NN effective interaction is rather weak (with 
a small increase with the increase of the energy and the 
angle) and, our conclusion is that it is not necessary to use 
it applying the folding approach to the cases of nucleon- 
nucleus scattering. 

iv) The results show that the spin-orbit interaction is 
rather important, particularly at angles larger than 60°, 
and that one can use in the Zs-term the microscopically 
calculated Re OP instead of the phenomenological WS 
potential with the three fitted parameters. 

v) The study of the dependence of the effective M3Y 
NN forces on the nuclear matter density shows small dif- 
ferences between optical potentials calculated with and 
without inclusion of the in- medium effect. The difference 
between the corresponding cross sections appears at larger 
angles and increases with the energy increase. 

vi) We showed that the effect of the non-linearity on 
calculations of Re OP connected with the factor jo (k(r)-s), 
where k(r) is the local momentum of motion, is small for 
the energy of 25.2 MeV/N but it increases with the energy 
increase up to an order of magnitude for i? c . m . > 40° for 
the energy E = 11 MeV/N. Thus, the non-linearity in the 
calculations of the Re OP should be taken into account in 
the calculations. 

Concluding, we would like to note that it follows from 
our results that the OP's can be calculated in the form 
Uopt = NrV f + iNjW H (i.e. with microscopically calcu- 
lated folding real part (V F ) and with calculated within 
HEA imaginary part (W )) using only two free param- 
eters (Nr and Nj) which renormalize the depths of the 
real and imaginary parts of OP. Thus, it is not necessary 
(at least on the basis of the existing experimental data 
on 6 He+p elastic scattering cross sections) to introduce 
a large number of fitting parameters, as is usually done 
in the arbitrarily chosen forms of the phenomenological 
and semi-microscopic optical potentials. It was pointed 
out that our approach can be applied to cases with en- 
ergies smaller than 100 MeV/N, like those of 71 MeV/N 
and 41.6 MeV/N considered in the present work. It can be 
concluded that microscopical optical potentials (includ- 
ing their imaginary part, e.g. of the HE A- type) rather 
than phenomenological ones have to be applied also for 
the cases of tens of MeV/N. Thus the approach can be 
used along with other more sophisticated methods like 
that from the microscopic g- matrix description of the com- 
plex proton optical potential and others. 
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